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ABSTRACT 

Composite laminates used in space structures are exposed to both thermal and mechanical loads. 
Cracks in the matrix form, changing the laminate thermoelastic properties. An analytical 
methodology is developed to predict microcrack density in a general laminate exposed to an 
arbitrary thermomechanical load history. The analysis uses a shear lag stress solution in 
conjunction with an energy-based cracking criterion. Experimental investigation was used to 
verify the analysis. Correlation between analysis and experiment is generally excellent. The 
analysis does not capture machining-induced cracking, or observed delayed crack initiation in a 
few ply groups, but these errors do not prevent the model from being a useful preliminary design 
tool. 


BACKGROUND 

The attractive properties of advanced composite materials make them ideal candidates for space 
structures. However, mechanical loads and wide swings in temperature cause cracks in the matrix, 
called microcracks. They can have a profound effect on the thermoelastic properties of the 
laminate, with potentially disastrous consequences in dimensionally critical applications. 

Prediction of microcracking under mechanical loading has been studied extensively. Many 
analyses have used in situ transverse ply strength as the cracking criterion. Lee & Daniel [1], for 
example, combined this failure criterion with a modified shear lag stress solution. Allen et al. [2] 
calculated damage with the internal state variable concept. Peters etal. [3] used a shear lag stress 
solution and Weibull strength distributions. However, Flaggs & Kural [4] questioned the validity 
of the strength-based approaches in general cases, showing that in situ ply strength is laminate 
dependent. Others have used fracture mechanics based cracking criteria. Laws & Dvorak [5] 
proposed a progressive damage model which incorporates a shear lag stress solution. Naim et al. 
[6] and Varna & Berglund [7] used variational solutions of the stress fields. Wang & Crossman 
[8] and Binienda et al. [9] calculated stress distributions with finite element models. Nuismer & 
Tan [10] combined a two-dimensional elasticity model with a fracture mechanics cracking criterion. 

Although many of the mechanical loading analyses have incorporated a residual thermal stress, 
very few predictive methodologies exist for progressive thermal loading. Most thermal analyses 
have focused only on the effects of microcracking on laminate properties. Bowles [1 1] and Adams 
& Herakovich [12], for instance, used finite element models to determine how coefficient of 
thermal expansion (CTE) is affected. Tompkins et al. [13] and Camahort et al. [14] 
experimentally determined the effects of thermal cycling on laminate thermoelastic properties. 
Some studies have been done to predict the onset of cracking (e.g. [15]). McManus et al. [16] 
developed an analytical method for thermally loaded cross-ply laminates. Crack density and 
reduced laminate properties were predicted as functions of monotonically decreasing temperature. 
The effects of thermal cycling were included in the analysis using a material degradation fatigue 
model. Experimentally, specimens were monotonically cooled or thermally cycled. Cracks were 
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subsequently counted on the ed~es. The analyses correlated very well with experiments. The 
work was later extended by P :< [17] to predict crack density in all plies of general angle 
laminates. Analytical predictK > agreed well with average interior crack densities measured 
experimentally. Park showed that, due to edge effects, interior crack counts are a better indicator 
of the laminate damage state than edge counts for other than cross-ply laminates. 


STATEMENT OF PROBLEM 

The objective of the present work is to develop an analytical methodology to predict the initiation 
and extent of microcracking in composite laminates under both thermal and mechanical loading. 
Analytical modeling in conjunction with experimental investigation is used to achieve this 
objective. An analysis is developed which, given material properties, laminate geometry, and load 
history, predicts crack densities and degraded laminate properties. Experimental work (i) provides 
verification of the analysis and (ii) allows a greater qualitative understanding of the problem. 
Finally, correlations provide insight into the details of the mechanisms by which microcracks form. 


ANALYTICAL METHODS 

Fig. 1 shows a laminate aligned with a global coordinate system xyz. The laminate is made up of 
unidirectional plies. Stacked plies with the same ply angle are assumed to act as a single thick ply, 
referred to as a ply group or layer. Cracks are assumed to span the layer thickness and propagate 
parallel to the fibers through the width of the laminate. Fig. 1 shows cracks of this type in a 
laminate with three ply groups. A local coordinate system x'y'z' is defined for each crack. The 
y‘ axis is aligned with the crack, parallel to the fiber direction of the layer, the x axis is aligned 
with the transverse direction of the layer, and the origin is at the center of the crack. 

To predict cracking in any one ply group, the laminate is modeled as being made up of two 
components: the cracking layer and the smeared properties of the rest of the laminate, henceforth 
represented by the subscripts c and r, respectively. A one-dimensional shear lag model, as shown 
in Fig. 2, is used to determine the stress and displacement distributions in the vicinity of a crack in 
its local x’y'z' coordinate system. In the shear lag solution that follows, E, a, u, and a are the 
stiffness, CTE, displacement, and normal stress, respectively, in the x direction. E„ is the 
effective laminate stiffness in the x direction. The shear stress between uncracked and cracked 
layers in the x'y' plane is q. The thicknesses of the cracking layer and the rest of the laminate are a c 
and a„ respectively, and a„ is the total laminate thickness. The laminate is subjected to a thermal 
load AT, which is the difference between the current temperature and some stress-free 
temperature, and an applied stress <r a , where 

a, = — (1) 


N x - is the laminate load in the x' direction. From laminate equilibrium. 


o a a 0 = <3,a r +o c a e 


From equilibrium of the cracked layer and the rest of the laminate, 

, = a .-£42i 

1 Alt " 


’ 2 dx' 

The stress-strain relations for the layers are 
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Substituting q = K(u e - u r ) (where K is an effective shear stiffness) into Eqns (3) and taking the 
derivative in x'. 
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Figure 1. Laminate geometry 

Combining Eqns. (4), multiplying by K, and substituting into Eqn. (5a), 

^ = f\f-T +{a ‘- a ' )Ar 
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Solving Eqn. (2) for a,, substituting into Eqn. (6), and rearranging. 
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where t, is the shear lag parameter, a dimensionless quantity that includes both material properties 
and geometry, and X. is used for convenience. From the rule of mixtures, 

a 0 E 0 = a r E r + a c E c (10) 

Substituting Eqns. (8), (9), and (10) into Eqn. (7) gives 
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This has a solution of the form. 
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Applying boundary conditions o c = 0 at x'= ±h gives 


A = 0 


4; 2 cosh(2 tyi/a e ) 

Placing Eqns. (13) into Eqn. (12), substituting Eqn. (9) into the result, and rearranging, 
_KaA (a r E r + a c E c )<y a , w-X cosh(2Cxya e ) 
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Solving Eqn. (2) for a,, substituting into Eqn. (14), and rearranging gives 

_<s.g, Ka‘ [“ fag, + a,g,)q. lArTl-^SfM 

' a, 2;’a,[ a.E,e, ( ‘ J { cosh(2?*M). 

r ' Substituting Eqn. (14) into Eqn. (4a), solving for d u c and integrating from 0 to x\ 


( 16 ) 


u - Ka cX' l (a r E r + a c E c )c a T g e sinh(2CsV a e ) 
‘ a,E r E 0 ' 2Cx’cosh(2 Qija e ) 


aATx' 


■ 2(?E c ^ E ' Ka { 1 2i'cosh(2^)JJ 

Substituting Eqn. (15) into Eqn. (4b), solving for du T and integrating from 0 to jc,' 
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The constants of integration in Eqns. (16) and (17) drop out because u,(0) = u £ (0) = 0. 

A new crack forms when the energy released due to crack formation, the strain energy release rate 
AG, is greater than the energy required to form a new surface, the critical strain energy release rate 
G, c . Assuming an existing uniform crack spacing 2 h, a new crack will instantaneously appear 
midway between two existing cracks when this failure criterion is met, as shown in Fig. 3. The 
minimum crack density at this load increment is defined as p, where p = \/{2h). 

G k is assumed to be a material property. AG is calculated from a Griffith energy balance, 

. — , AW -ALT 

AG (18) 

a c 

where AW and AU are the changes in external work and internal strain energy, respectively, 
between the states before and after the hypothetical new crack forms. The strain energy 
contributions from normal stresses, and shear stresses, £/,, are 


u °=\& 




The total strain energy is U = U a +U r The change in strain energy when a new crack appears is 
found by subtracting the strain energy of the representative volume before the crack forms, where 
the cracks are separated by 2 h, from that after formation, where the cracks are h apart: 


tM =[2U\ k -U\»} 


Solving Eqn. (20) yields 


AU = ~f r a ‘ E 'vi at ~ a ^ A7 ~ 2 2tanhf —1 - tanhf (21) 

Note that there is no coupling between the thermal and mechanical loading in Eqn. (21). However, 
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Figure 2. Shear lag model 


Figure 3. Geometry used 
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the work done by the applied loading is 


W = a„a 0 u r 

The change in external work when a crack forms is given by 

&W = [2W] t -W\, k ] 

Placing Eqn. (17) into Eqn. (22) and solving Eqn. (23), 
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Placing Eqns. (21) and (24) into Eqn. (18) gives the total change in strain energy release rate: 


AG = 
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It is important to note that the final expression for AG has a thermomechanical interaction term due 
to the coupling between thermal and mechanical loading in Eqn. (24). Note also that Eqn. (25) 
must be solved numerically or graphically to find h (or p) for a given AT and cr a . 


At every load increment, the crack formation analysis is carried out for each ply group in turn. The 
properties of the ply group considered are calculated in its local coordinate system x'y'z' using 
familiar CLPT relations, and knockdown factors to account for existing cracks. Next, Eqn. (25) is 
used to determine the new crack density in that layer. Once all of the ply groups have been 
examined at this load increment, the effective laminate properties are updated. The cycle is 
repeated for each load increment until the entire load history is completed. This iterative damage 
progression model is implemented by the computer code CRACKOMATIC H, available by request 
from the authors, which takes material properties, laminate geometry, and thermomechanical load 
history and outputs crack density and degraded laminate properties. The analysis and its 
implementation are presented in detail in ( 1 8]. 


EXPERIMENTAL PROCEDURES 

Experiments were conducted to correlate with analytical predictions and to investigate the 
microcracking problem qualitatively. AS4/350I-6 graphite epoxy laminates were fabricated at the 
Technology Laboratory for Advanced Composite Materials at MIT using 0.127 mm (0.005") thick 
plies. Three different layups, [04/454/904/-454]s, [0 j/ 452/902/-45 2 ] s , and [(y602/-60 2 ] s , were 
manufactured to investigate laminate geometry and thickness effects. Panels were cut into test 
specimens according to the matrix in Table 1. Both long edges of every specimen were polished to 
minimize crack initiation sites and to facilitate microscopic inspection. All specimens were dried 
and then stored in desiccant to eliminate moisture. 

Thermal specimens were monotonically cooled to progressively lower temperatures in an 
environmental chamber. Starting at room temperature, die specimens were cooled to temperatures 
as low as -184°C (-300°F) at a rate of 14°C/min. The cooling and heating rates were low enough to 
avoid significant thermal gradients. Specimens were soaked at the target temperature for 5 minutes 
and returned to room temperature. The edges of the specimens were then inspected for cracks 
under an optical microscope at lOOx magnification. Sacrificial thermal specimens were sanded to 


Table 1. Test matrix repeated for all 3 layups 


Tesring/Daia 

Dimensions 

(cm) 


Thermal/Edge 

7.62 x 2.54 

5 

7.62 x 1.27 

5 

Thennal/Interior 

7.62 x 2.54 

15 

7.62 x 1.27 

15 

Mechanical/Edge 

25.4 x 2.54 

5 


Table 2. Material properties 


Ei (GP» 

142.0 

tkm 2E3B 

-0.36 

E, (GPa) 

9.81 


28.8 

V 

0.30 

EF8EE Sil 

241.0 

G„ (GPa) 

6.0 

c 

1.0 



























incremental depths to measure interior crack densities. An observed crack extending more than 
half the thickness of a ply or ply group was counted as a crack. Mechanical coupons were loaded 
under monotonically increasing tension at room temperature. Edge crack data was collected at 
incremental loads until failure or other damage modes, such as delamination, were observed. 


RESULTS 

Material data used in the analyses are presented in Table 2. The same values of the shear lag 
parameter, £, and the critical strain energy release rate, G te , were used for all analyses. In each of 
the following figures, the measured and predicted crack densities for specified ply groups are 
shown as functions of the mechanical or thermal load. The results for the 45« and 90 4 ply groups 
in mechanically loaded [0 4 /454/90 4 /'45 4 ] s laminates are shown in Fig. 4. The 45* ply group begins 
cracking at a higher load and accumulates less cracks than the 90« layer. The scatter is relatively 
low, except in the region where cracking initiates. The analysis successfully predicts the onset of 
cracking and only slightly underpredicts the accumulation of cracks in both layers. The same 
trends are seen in the 45 2 and 90 2 ply groups in the [O2/452/9O2/-45JS laminates, shown in Fig. 5. 
Note that cracking initiates at nearly identical loads in both laminates, even though the latter are half 
as thick. Additionally, the thinner [(y452/902/-45 2 ]s laminates reach significantly higher crack 
densities than the [04/454/904/- 4 5 4 ] s laminates. 

The analysis does not predict the initiation of cracking in the middle ply groups under mechanical 
loading satisfactorily. Cracking initiates in these layers only when the adjacent layers start to 
crack, even though the analysis predicts cracking much earlier. However, the analysis appears to 
capture microcrack accumulation once cracking has begun, as shown by the [(VftV-tiOJj results 
in Fig. 6. Here cracking is suppressed in the -60 4 layer until the load reaches 13kN, when 
cracking initiates in the adjacent 60 2 layers. Once the -60« layer begins to crack, the analysis 
predicts crack density in that layer nearly perfectly. Correlation between analysis and experiment 
in the 60? layer is very good throughout the entire range of microcracking. 

Thermal loading data and analyses for the 45 4 and 90 4 ply groups in the [0 4 /45 4 /90 4 /-45 4 ] s 
laminates are shown in Fig. 7. Scatter in the thermal data is generally higher than in the mechanical 
data. The analysis captures the cracking behavior of the 45 4 layer, including initiation temperature 
and crack accumulation. Some initial cracks were present in the 90 4 ply group after manufacture, 
which the analysis does not predict. Nevertheless, the analysis appears to successfully predict 
crack density in this layer at the lowest temperature increments. 

The thermal data does not manifest the middle ply problem described earlier for the mechanical 
data. Cracking initiation in the middle layers is not influenced by adjacent layers. For example, 
the analysis predicts initiation and accumulation in the -45j ply groups in the [0 4 /45 4 /90 4 /-45 4 ] s 
laminates very well, shown in Fig. 8. Fig. 9 shows the thermal loading data and analyses for the 
45 2 and 9O2 ply groups in the [0^45^0^-45^3 laminates. Little cracking was observed in these 
laminates under thermal loading. The analysis predicts no cracking down to -184°C, agreeing well 
with experiment, as virtually no cracks were observed in the 45 2 layer. The onset of cracking in 
the 90 2 layer is later than predicted, but the analysis agrees well with the data at the lowest 
temperature. Unfortunately, it is difficult to assess the correlation between experiment and analysis 
for this layup, as so few data points are available. The [0 2 /60y'-60 2 ] s laminates do not crack at all 
under thermal loading, even though the analysis predicts initiation in all layers. 


DISCUSSION AND CONCLUSIONS 

1. The present methodology shows promise as a powerful, unified analytical tool for predicting 
microcracking in all plies of a general angle laminate exposed to an arbitrary thermomechanical 
load history. It could be useful in any application of composites where dimensional stability is 
critical and the structure is exposed to a combination of thermal and mechanical loads. 

2 . Its general applicability is further bolstered by the fact that the shear lag and fracture mechanics 
parameters do not demonstrate laminate dependence to an acceptable approximation. Thus they 
only need to be measured once for a given material system. This can be done by counting 
cracks in a cross-ply laminate under tension and using the analysis to back out the parameters. 


3 . The limits of the analysis need further exploration. Different material systems, or very thick or 
very thin ply groups, may have different mechanisms of damage that fall outside of the 
assumptions used in the model presented here. 

4. The analysis assumes that critical starter cracks are inherent in the uncracked laminate. If these 
flaws are edge defects, then laminates with a strong edge effect (see [18]) may experience 
delayed onset of cracking. Cracking in a layer is then initiated instead by stress concentrations 
from cracks in adjacent layers. This offers an explanation for the delayed onset of cracking in 
some layers under mechanical loading. This could also explain why the [O 2 / 6 O 2 /- 6 O 2 JS 
laminates, which have a strong edge effect, did not crack under thermal loading. 

5. Even though the analysis does not capture crack initiation in some layers under mechanical 

loading, it predicts crack accumulation very well after the onset of cracking. The analysis also 
sometimes predicts cracking when none are observed experimentally, such as the thermal 
loading of the laminates. The methodology is conservative in these cases, and 

works well in all others considered. Hence it appears to be a very useful analytical tool. 
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Figure 5. Mechanical loading of 16-ply 
specimens 
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Figure 6. Delayed initiation of cracking in Figure 7. Thermal loading data in 32-ply 
middle layer under mechanical loading specimens 
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Figure 8. Thermal data; 32-ply specimens Figure 9. Thermal data; 16-ply specimens 








